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ABSTRACT 

In this article the comparison of the dynamic response of a cam follower 
system is made between several laws of motion. The laws that are compared are 
the study of a hypothesis that is intended to valídate denominated Bezier Degree 
15, against standard laws used in the design of cams like polinomial law 3-4-5, 
polinomial law 4-5-6-7, the cycloidal law and the Bezier in stages. For this purpose 
three different dynamic models are used: a model of a degree of freedom that 
supports the residual vibration, the model that considers the flywheel to verify 
torsional vibration and a model of two degrees of freedom to identify the takeoff of 
the follower. Once the comparative advantages are demonstrated, the cam design 
and its subsequent construction are carried out. Then a study of the generated 
noise is made to support the proposal presented. 
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APLICACIÓN DE CURVAS DE BÉZIER EN EL DISEÑO 
Y OPTIMIZACIÓN DE LEVAS PARA ALTA VELOCIDAD 


RESUMEN 

En este artículo se efectúa la comparativa de la respuesta dinámica de un 
sistema leva seguidor entre varias leyes de movimiento. Las leyes que se 
comparan son la hipótesis de estudio que se pretende validar denominada Bézier 
Grado 15, frente a leyes estándar utilizadas en el diseño de levas como la ley 
polinomial 3-4-5, la ley polinomial 4-5-6-7, la ley cicloidal y la Bézier a tramos. Se 
recurre para tal fin tres modelos dinámicos diferentes que son: un modelo de un 
grado de libertad que sustenta la vibración residual, el modelo que considera el 
volante de inercia para verificar vibración torsional y un modelo de dos grados de 
libertad para identificar el despegue del seguidor. Una vez evidenciadas las 
ventajas comparativas se efectúa el diseño de la leva y su posterior construcción. 
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Luego se realiza un estudio del ruido generado para mayor sustento a la 
propuesta presentada. 

Palabras clave: Bézier, Levas, Modelamiento dinámico. 


1. INTRODUCCIÓN 

Actualmente el laboratorio de mecanismos de la Universidad de las 
Fuerzas Armadas cuenta con el equipo de análisis de levas Tecquipment TM 21, 
el cual data del año 1979. Este equipo ha sido rehabilitado y automatizado, en 
trabajos de investigación previos en el aspecto de la adquisición de datos [1]. Sin 
embargo, la parte mecánica, es decir las levas originales que se utilizan en dicho 
equipo no han sido renovadas, creándose un desfase entre la teoría y la práctica 
de laboratorio debido a que estas levas eran diseñadas mediante una sucesión 
de arcos tangentes entre sí y probablemente fabricadas con maquinaria 
convencional, lo que conlleva a que las curvas de aceleración del seguidor sean 
siempre discontinuas, ver Fig. 1 [2], 



Figura 1 . Desplazamiento, velocidad y aceleración de las levas del TM21 

Este hecho viola el principio del diseño de las levas [3] que se evidencia 
por el ruido del golpeteo generado. La teoría matemática de las levas ha 
evolucionado de la mano de los procesos de manufactura pudiéndose encontrar 
en la literatura especializada un sinnúmero de leyes de desplazamiento, las 
cuales son estructuradas a tramos para generar el perfil de la leva, lo que 
implicará que alguna de las derivadas del desplazamiento sean discontinuas. La 
hipótesis que se plantea en el presente trabajo es utilizar para el perfil, las curvas 
de Bézier. La curva de Bézier fue desarrollada por Pierre Étienne Bézier para 
reconstruir matemáticamente las formas libres de las carrocerías y tiene la 
característica de ser globalmente deformable y predecible [4]. La curva de Bézier 
será aplicada en este estudio para garantizar que la ley del seguidor subida- 
reposo-bajada y todas las funciones cinemáticas, velocidad, aceleración, sobre 
aceleración y demás derivadas sean continuas, con lo que se tendrá un 
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funcionamiento suave y silenciosos. Esto va a ser validado mediante una 
comparativa entre la hipótesis y las leyes usuales con el objeto de construir la 
leva mencionada. En el estudio comparativo se utilizan tres modelos dinámicos y 
se resuelven tres ecuaciones diferenciales, ya que cada una de ellas toma en 
consideración diferentes parámetros, que definen vibración, vibración torsional y 
despegue del seguidor. Para esto se utilizará el programa de análisis 
computacional Wolfram Mathematica. Finalmente se construirá la leva con el 
perfil de Bézier, con el fin de comprobar de manera práctica los resultados 
obtenidos en el estudio dinámico. Se pretende finalmente que los estudiantes que 
toman el curso de mecanismos perciban el vínculo entre teoría y práctica en lo 
referente a la práctica de levas. 

2. GENERACIÓN DEL PERFIL 

En el presente trabajo se realizó un estudio comparativo entre la hipótesis 
propuesta que se denominará Bézier Grado 15 y cuatro leyes de desplazamiento 
estandarizadas: polinomial 3-4-5, polinomial 4-5-6-7, cicloidal y Bézier a tramos 
[5], 

Una curva de Bézier paramétrica de grado n es una combinación lineal de 
polinomios de la base de Bernstein BJ 1 (u) sobre un dominio unitario, Ec. (1) 

Dónde: 

n= grado del polinomio; ¿= contador de 0 a n 

Esta operación genera un vector, cuya transpuesta a su vez es 
multiplicado por un polígono de control P. Para generar la curva de Bézier grado 
15 se utilizó la herramienta iterativa Manipulate[ ] del software Wolfram 
Mathematica, con la cual se sintonizo el vector puntos de control al requerimiento 
subida-reposo-bajada. Por tanto, el polígono de control que se utiliza es: 

P T = [0,0,0,0,1.5,15,0.868,0.868,0.868,0.868,1.5,1.5,0,0,0,0]. 
Obteniéndose la curva Bézier Grado 15 de la Ec. (2) correspondiente a la Fig. 2. 

m = -2047.5 (?) 4 (f - l)“ + 2047.5 (!)” (? - l)‘ - 4344.34 - l)’ + 

4344.34 (£)’ (í - 1)‘ 4 4504.5 $ - l)“ - 4504.5 (J)“ (í - 4)' 4 

5585.58 (í) 7 (í-l)'-5585.58 (0“(í-l) 7 

( 2 ) 
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Figura 2. Ley Bézier Grado 15 propuesta 

Cabe aclarar que esta ecuación representa subida-reposo alto-bajada, 
descartándose incluir el reposo bajo debido a la complejidad de la ecuación 
resultante, complejidad que no tendría sentido si se busca una utilidad práctica del 
diseño propuesto. La ecuación final por tanto es la Ec. (3). 

= fs(0)> 0<6<p 

K ’ 1 0 , p < e < 360 ° 


3. MODELAMIENTO DINÁMICO 


• Vibración residual 

El modelo matemático de un grado de libertad se basa en la configuración 
del equipo de levas TM 21 del laboratorio de mecanismos de la Universidad de las 
Fuerzas Armadas como se observa en la Fig.3. 


d<*J ]_1 k1{x) 
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Figura 3. Equipo TM21 y sistema masa resorte 


Se supone que k2 » kl y que x > y 
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^ F = m x 

= — /qx — c x x— k 2 (x — y) = mx (4) 

La ecuación diferencial a resolverse es por tanto Ec. (5): 

mx+ c t x + (k 1 + k 2 )x = k 2 y (5) 

Dónde: 

m = masa del seguidor; ci = amortiguación hidráulica; ki = rigidez del resorte; = 
rigidez del seguidor; y = ley del seguidor 

Para resolver la Ec. (5) es necesario introducir los parámetros físicos e 
inerciales definidos previamente, Tabla I 


Tabla 1: Parámetros inerciales [5] 


Nombre de la propiedad o parámetro 

Valor 

Constante de resorte real, kl en [N/m] 

4030 

Rigidez del actuador/seguidor, k2 en [N/m] 

330000000 

Amortiguamiento, c en [N*s/m] 

8.642 

Masa del sistema, m [kg] 

1.3 

Altura/alzada de la leva, h [m] 

0.05 


Como también las ecuaciones del desplazamiento de las diferentes leyes 
de desplazamiento utilizadas en este estudio. Se integran las ecuaciones 
mediante el software Wolfram Mathematica utilizando el comando NDSolve[ ] que 
resuelve numéricamente la ecuación diferencial, Evaluatef ] para evaluar el 
resultado y Plot[] para visualizar gráficamente los mismos con la cual se puede 
realizar la comparación. 

En la Fig. 4 se puede observar los resultados de la solución de la Ec.5, 
donde se compara el comportamiento dinámico de las cinco leyes del seguidor 
utilizadas. 



— Bézier Grado 15 

— Polinomial 345 

— Polinomial 4567 

— Cicloidal 

— Bézier a tramos 
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Figura 4. Resultado gráfico de desplazamientos EDO 1 grado libertad 


Para observar de manera más clara y objetiva cuál de los perfiles tiene 
mejor comportamiento dinámico o menor vibración residual se genera la Fig. 5, la 
cual es un acercamiento del tramo correspondiente al detenimiento en la bajada 
del seguidor 



Bézier Grado 15 

— Polinomial 345 

— Polinomial 4567 
— Cicloidal 

— Bézier a tramos 


Figura 5. Vibración residual en el reposo bajo 
La curva de aceleración de la Fig. 6 es directamente proporcional a las 
fuerzas de inercia que percibe el seguidor, claramente se observa que la curva 
más ventajosa es la propuesta en el presente estudio Ec.(2) tanto por la menor 
cantidad de oscilaciones como por su menor valor relativo. 



• Vibración torsional 

El equipo TM21 posee un volante de inercia que provee estabilidad al 
sistema, controlando las fluctuaciones de velocidad, ver Fig.7, en vista de que por 
lo tanto se inducen vibraciones torsionales se efectúa en este apartado la 
correspondiente comparativa 
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Figura 7. Esquema del equipo de levas con volante de inercia 

Así pues, las ecuaciones de movimiento se pueden obtener de la 
sumatoria de pares de torsión [6]. 

2 T = I e ( 6 ) 

Que desarrollados da, Ec.(7): 

ie c + k s [e c (!)-e l (t)] = T leva (7) 

Dónde: 

I = momento de inercia de masa; 0 C = aceleración angular de la leva; ks = rigidez 
del eje; 0 C = posición angular de la leva; 0¡ = posición angular del volante de 
inercia; T leva = Par torsional de la leva 

Para obtener el Torque externo se iguala la potencia de entrada con la 
potencia de salida [7] en el sistema leva seguidor, Ec. (8). 

Tleva Óc = F cX ( 8 ) 

Dónde: 

F c = Fuerza de contacto; x = velocidad del seguidor; 0 C = velocidad angular de la 
leva 


La fuerza de contacto resulta de la sumatoria de fuerzas actuantes sobre 
la leva, ver Fig. 8. 



Figura 8. Esquema de la leva y diagrama de cuerpo libre 
La fuerza cinetostática de contacto está dada por Ec.9: 
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F c (t) = —m x — tic — kfX (9) 

Reemplazando la Ec. (9) en la (8) y (7) se obtiene la siguiente ecuación 
diferencial, Ec. (10) 

X 

I é c + k s [0 c (t) — + (m x + ex + kfx)— = 0 (10) 

e c 

Los valores que se grafican en la Fig. 9 corresponden a la diferencia que 
existe entre la posición angular de la leva y la posición angular del volante de 
Inercia: 



Bézier Grado 15 

— Polinomial 345 

— Polinomial 4567 

— Cicloidal 

— Bézier a tramos 


Figura 9. Vibración torsional 

La gráfica muestra que la leva Bézier Grado 15 exhibe mejor control de la 
vibración torsional. 

• Modelo dinámico para predecir el despegue del seguidor 

Es importante añadir el análisis de despegue entre el seguidor y la leva, ya 
que de esta manera se determina la ley que mejor resiste el incremento de 
velocidad angular en la leva sin evidenciar despegue del seguidor. El modelo que 
predice el salto o despegue, es el modelo de dos grados de libertad según la Fig. 
10 . 



Figura 10. Modelo de dos masas simplificado 

Este modelo matemático toma en cuenta que la masa del sistema está 
dividida entre el seguidor y el efector final. Para plantear la ecuación diferencial se 
puede usar el método de ensamblaje directo usado en el estudio de elementos 
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finitos, donde se considera mi, m2, m3 las masas en los nodos del sistema de 
resortes por tanto las matrices de rigideces y amortiguación parciales son: 

Nodos 1 y 2 Nodos 2 y 3 



1 


Cl 

-q 




Efectuando el ensamblaje se obtiene, Ec. (11): 


' kl 

~h 

o ■ 

-x x 


' Ci 

-Ci 

0 ■ 

-Xi 


m 1 

0 

0 

X! 


\ F c] 

-fcl 

k x + k 2 

-k 2 

*2 

+ 

-Ci 

Ci+C 2 

~ c 2 

*2 

+ 

0 

m 2 

0 

x 2 

= 

0 

. 0 

~^2 

^2 - 

* 3 . 


. 0 

~ c 2 

c 2 - 

-X 3- 


. 0 

0 

m 3 . 

-x 3 - 


.0. 


Dónde: 

ki = rigidez del vástago del seguidor; = rigidez del resorte de cierre de fuerza; 
ci = amortiguación del vástago del seguidor; C 2 = amortiguación viscosa del 
seguidor con respecto a tierra; mi = masa equivalente concentrada en el lado del 
rodillo; rri 2 = masa equivalente concentrada en el lado del efector; rri 3 = masa 
ficticia de la referencia fija; F c = fuerzas de contacto; s = ley de la leva; z, y = 
respuestas dinámicas de mi y nri 2 respectivamente 


Es claro que ki» k 2 

Considerando que el nodo correspondiente a la masa m3 es fijo, el 
sistema de ecuaciones diferenciales queda, Ec. (12): 


fci -k x 1 rXij r Ci -c x 1 rxjl fm 1 0 1 pé'il = n 

.—k-i k-i + J ^2-1 L ^1 Cl ~t~ ^2-1 lx 2 J - 0 77*2-1 \-%2- *- 


( 12 ) 


Usando los parámetros y y z de la Fig. 10 y tomando en cuenta que la 
fuerza externa en la masa mi es la fuerza de contacto F c (t), Ec. (13). 


' k x —k x 
-k x k 1 + k 2 


r z i + r Cl ~ Cl i[ ¿ l. K 

LyJ 1-C X c x + c 2 i [y\ + o 


Finalmente desarrollando la Ec. 
diferenciales Ec.(14) y (15). 


o 

m 2 


( 12 ) 



se obtienen dos ecuaciones 


c x ¿ — c t y — k t y + k t z + m^z = F c (14) 

y{c 1 + c 2 ) + z(k 1 + fe 2 ) — c x z — k 2 z + m 2 y = 0 (15) 

Las cuales se resuelven de la siguiente manera: 

Inicialmente se considera que el seguidor y la leva no se separan, por 
tanto z = s, siendo s la ley del seguidor, por tanto, se resuelve numéricamente la 
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Ec.15, reemplazando la z por s como si se tratara de un sistema de un grado de 
libertad, Ec. 16, con condiciones iniciales cero. 

y(c ± + c 2 ) + s(k 1 + k 2 ) — c x s — k 2 s + m 2 y = 0 (16) 

De la solución de Ec. (16) se obtiene y, y. Luego se evalúa para cada 
tiempo la fuerza de contacto (Ec. 15) según el formato indicado en la Ec.17 


c-^s — — k t y + k-tS + m^s = F c (17) 

Se determina el instante ti en que esta fuerza se hace cero, lo que implica 
el inicio del despegue. De la solución de Ec. 16 se determina el desplazamiento 
y(tl) y la velocidad y (t ± ), también se considera que z (t x ) = s (t x ) y z (t x ) = 
¿(tj), siendo estas las cuatro condiciones iniciales necesarias para resolver el 
sistema homogéneo de dos grados de libertad, Ec. 18 y 19: 

c !¿ — c iy — k t y + k t z + m-^z = 0 (18) 

y(c ± + c 2 ) + z(k ± + fc 2 ) — c x z — k 2 z + m 2 y = 0 (19) 

En esta segunda solución se testea z contra s, cuando z < s , nuevamente 
se establece el contacto y nuevamente se debe usar la Ec. (16) con la (17). Una 
vez resuelto el modelo se puede observar la gráfica de la respuesta del sistema de 
ecuaciones diferenciales, Fig. (11), aquí se puede identificar una ligera 
prominencia al momento de finalizar la etapa de subida. Esta prominencia es la 
evidencia del despegue del seguidor. 


BeziErGrada 15 


Palinomial 345 


Bezier a tramas 






Palirmmial 45G7 


Cicloidal 






Figura 11. Respuesta dinámica del despegue 

Finalmente se determinó la sensibilidad al despegue en función de la 
velocidad angular para los cinco perfiles comparados, Tabla II. 


Tabla 2: Sensibilidad al despegue 
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Leva 

Velocidad 

[rpm] 

Tiempo 

[s] 

Bézier Grado 15 

570 

0.9 

Polinomio 345 

415 

1.5 

Polinomio 4567 

390 

1.38 

Cicloidal 

410 

1.38 

Bézier a tramos 

400 

1.62 


Se constata nuevamente que el perfil Bézier Grado 15 es la que más 
retarda el despegue del seguidor con respecto a la velocidad. 

4. FABRICACIÓN DE LA LEVA PROPUESTA 

Una vez comprobadas las características satisfactorias del perfil propuesto 
se procede a desarrollar el modelo CAD, utilizando el software Mathematica para 
generar el perfil de la leva según la Fig. 12. 



Este perfil posteriormente será convertido a una región y luego a un archivo *.dxf. 
En la Fig. 13 se observa el modelo final y el Postprocesado CAM 
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Para su mecanizado final según Fig. 14. 



Figura 14. Leva final 


5. PRUEBAS DE FUNCIONAMIENTO 

• Comparaciones previas 

La prueba de funcionamiento de la leva se la realizó con el fin de constatar 
el desempeño de la misma, en lo que al ruido generado se refiere, al ser sometida 
a altas velocidades en el equipo de levas del laboratorio de mecanismos 
Tecquipment TM21. 

Para lo cual se utilizó un sonómetro digital JTS1357 con capacidad de 30 
a 130 dB, ver Fig. 15. 
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Figura 15. Sonómetro digital 

Generando los siguientes resultados, Tabla III. Donde se muestra los valores 
obtenidos de ruido de la leva Bézier Grado 15, la leva propia del equipo original 
TM21 y una leva Polinomial 345. 

Tabla 3: Tabulación del ruido a distinta velocidad 



Bézier Grado 15 

Leva Equipo TM 

21 

Polinomial 

345 


RPM 

dB 

RPM 

dB 

RPM 

dB 

40 

145 

73.4 

159 

70 

158 

72.5 

45 

180 

78.2 

188 

72 

187 

76.4 

50 

210 

79 

220 

74 

219 

80.4 

55 

240 

79.2 

250 

75.7 

248 

81.4 

60 

280 

81.3 

285 

77.6 

273 

82.4 

65 

310 

81.5 

313 

80.6 

308 

83.9 

70 

345 

82 

353 

82.8 

342 

85.3 

75 

395 

86.5 

388 

94.7 

384 

88.7 


De la Tabla III se efectúa la regresión polinomial lineal 


Bezier = 25.93808 + 0.57355 x - 0.002045 x 2 + 0.000002483 x 3 (20) 

Leva Equipo TM21 = 0.41488 + 0.84929x - 0.0033286 x 2 + 0.000004537x 3 (21) 
Polinomial = 2.293766 + 0.78001x - 0.002606 x 2 + 0.000003025 x 3 (22) 


Se observa en la Fig.16 que a partir de los 335 rpm la leva Bézier Grado 
15 propuesta tiene más nivel de ruido, antes de esta velocidad la leva cóncava 
original del equipo tiene menor nivel en virtud de que no posee reposo alto. 

Sin embargo, después de la velocidad indicada, la ley Bézier Grado 15 
muestra su ventaja en la reducción de ruido. 
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Figura 16 . Gráfica de valores obtenidos del sonómetro digital 


CONCLUSIONES 


El principal objetivo de este trabajo que fue dotar al equipo TM21 de una leva 
acorde a la teoría actual de Bézier específicamente desarrollando una curva 
paramétrica de grado 15 de tal manera que la teoría y la práctica de laboratorio 
guarden concordancia. Si bien la leva de Bézier ofrece mejor respuesta dinámica 
en todos los aspectos, su uso se vería limitado por el escaso control que presenta 
en la temporización exacta de los eventos. Se podría profundizar el proyecto para 
resolver la incertidumbre de que la mejora observada es causa de la continuidad 
de la ecuación y/o al hecho de que las subidas o bajadas presentan menos 
pendiente y por lo tanto las aceleraciones disminuyen drásticamente. Se verificó 
que la leva diseñada ofrece mejor control al ruido sobre los 335 rpm comparada 
con las levas construidas con las leyes convencionales. El presente trabajo fue 
desarrollado usando el método de los parámetros concentrados, se podría 
continuar la investigación usando el método de elementos finitos usando análisis 
estructural transitorio. También se evidenció que las curvas polinomiales de más 
alto grado tienen mejor control de vibración, pero son más sensibles al despegue 
del seguidor. 
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